Travelling to exotic places with cavity QED systems 
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Recent theoretical schemes for utilizing cavity QED models as quantum simulators are reviewed. 
By considering a quadrature representation for the fields, it is shown how Jahn- Teller models, 
effective Abelian or non-Abelian gauge potentials, transverse Hall currents, and relativistic effects 
naturally arise in these systems. Some of the analytical predictions are verified numerically using 
realistic experimental parameters taking into account for system losses. Thereby demonstrating 
their feasibility with current experimental setups. 
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I. INTRODUCTION 

Back in 2006 I was a co-author on a paper called "Trav- 
elling to exotic places with ultracold atoms " where we dis- 
cussed some, at the time, new ideas how to mimic unusual 
phenomena occurringin condensed matter systems by us- 
ing ultracold atoms [f| ■ Due to the controllability of sys- 
tem parameters and long decoherence times, cold atoms 
in optical lattices had shown to be suitable candidates for 
simulation of different many-body theories developed in 
the community of condensed matter physics [2j . Around 
the same time, many theoretical predictions were pre- 
sented in which high energy and cosmological phenom- 
ena could be simulated by means of Bose-Einstein con- 
densates [3J. In addition, ion-trap systems have recently 
drawn much attention in terms of simulating various spin 
chains as well as relativistic effects 

In the present proceeding I summarize some new direc- 
tions employing cavity QED systems as quantum simula- 
tors. On the experimental side, cavity QED systems have 
seen great progress during the last years. Pure quan- 
tum effects arc nowadays not restricted to single few-level 
atom-cavity setups [5[, but include both quantum-dots 
(qdots) or SQUfDs [6j as well as Bose-Einstein conden- 
sates [3] coupled to single cavity modes. These new ex- 
periments pave the way for reaching the superstrong cou- 
pling regime where the strength of the effective matter- 
field coupling exceeds the system losses by several or- 
ders of magnitude. One crucial outcome of reaching this 
regime is that application of the rotating wave approxi- 
mation (RWA) becomes questionable [8]. In most cases, 
going beyond the RWA implies that the corresponding 
system Hamiltonian is not integrable and other approx- 
imations or numerical methods are required. As out- 
lined in Ref. 0, by expressing the cavity field(s) in its 
quadrature operators, it directly follows that beyond the 
RWA typical cavity QED models describe a set of cou- 
pled harmonic potentials. In this representation it is 
thereby possible to develop new physical intuition not 
easily extracted when the field is expressed in boson cre- 



ation and annihilation operators. Moreover, in this pic- 
ture the link between cavity QED models and systems 
from other fields of physics readily follows. 

We consider a single iV-level atom |9j coupled to a set 
of quantized cavity modes. The quadrature operators for 
mode k reads 



i 

V2 



(a fe + 4), P k = ^ (a fc -4) (I) 



with dtk (fit) the annihilation (creation) operator of a 
photon in mode k. The quadrature operators obey the 
canonical commutation relations [Xfe,P;] = idki, and 
thereby can be viewed as "position" and "momentum" 
for the field. In the dipole approximation we assume the 
electric field -E(x) to be constant over the extent of the 
atom, hence letting x = 0, to obtain 



E — ^ £k£kPk, 



(2) 



where e~k and £k are the polarization vector and the field 
amplitude of the fc'th mode respectively. Denoting the 
dipole moment between the atomic states \i) and \j) by 
da, the atom- field interaction takes the form 



(3) 



Due to selection rules or large atom-field detunings, 
most terms in the above sums vanish either exactly or 
approximately. Explicitly, the dipole moment d%j = 
(<•</;',•'/,'.,:•• with df 3 = -e|i)(«|a|j)0'| +H.c and e the 
electron charge. The full Hamiltonian, containing the 
free field and internal atomic energies, becomes 



H — Hf + H a + Hi, 



(4) 



where 



H f = hj2^(^-+^), Ha^mui (5) 
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LOk the frequency of mode fc, and Ei being the energy of 
the j'th initial atomic state. 
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II. JAHN-TELLER SYSTEMS 

Jahn- Teller systems frequently occur in both molecu- 
lar/chemical physics and condensed matter physics [Iol |. 
The name dates back to the works by Hermann Jahn 
and Edward Teller studying symmetry breaking for poly- 
atomic molecules . Common for the Jahn- Teller mod- 
els is that two potential surfaces become degenerate in a 
single point forming a conical intersection. About two 
decades after the pioneering work by Jahn and Teller, 
Longuct-Higgins showed that by encircling the degener- 
acy, the electronic wave function changes sign [l2] . This 
phase factor was later realized to be the so called Berry 
phase |13[. 



A. The P x E model 



been known in cavity QED, and in the case of infinitely 
many atoms (the thermodynamical limit) it describes a 
second order quantum phase transition [15j |. However, 
only recently has it been shown that this phase transition 
is identified as a Jahn- Teller symmetry breaking [161 ]. 




The simplest non-trivial atom-cavity system considers 
a two-level atom interacting with a single cavity mode. 
Introducing the Pauli matrices a x — |1)(2| + |2)(1|, a y = 
i(|l)(2|-|2)(l|), and u z = |1)(1|-|2)(2|, the Hamiltonian 
reads 

( P 2 X 2 \ Ml 
Hps = h\^> I — + — J + — o- z + hga x P. (6) 

Here, Q is the atomic transition frequency and g the 
effective atom-field coupling. Application of the RWA 
with respect to the first two terms renders the well 
known Jaynes-Cummings model For true atomic 

systems, it accurately explains several experimental ob- 
servations @. In terminology of molecular physics, HpE 
is equivalent to the (3xE Jahn- Teller Hamiltonian by sim- 
ply interchanging X and P. For a deeper understanding 
of the dynamics, it is more convenient to think of X 2 /2 
as the kinetic energy term and P as the coordinate. In 
this picture, the adiabatic potentials are defined as 



hwP 2 



±h 



n 2 

1 r+9 2 P 2 - 

4 J 



(7) 



The diabatic potentials are V^(P) = huP 2 /2±hgP. Fig- 
ure [1] depicts examples of the different potentials. When- 
ever \g\ > the lower adiabatic potential V~JP) 
possesses two minima, while for \g\ < yuQ/2 a single 
global minimum occurs at P = 0. The Jahn- Teller ef- 
fect [HI EH is easily understood from the adiabatic po- 
tentials; the amplitude of the ground state wave function 
does not attain its maximum at P = 0, but around the 
two shifted minima. For polyatomic molecules it implies 
that the conical intersection cause the molecule to vibrate 
around a less symmetric atomic configuration, or in other 
words, the ground state energy is decreased by lowering 
the symmetry. In terms of cavity QED, the Jahn- Teller 
effect states that for sufficiently strong atom-field cou- 
pling the ground state is not the one with the atom in its 
ground state and the field in vacuum. This effect has long 



FIG. 1: Typical diabatic (dashed) and adiabatic (solid) po- 
tentials for the f3 x E model. 

Before proceeding, some words of caution. The pres- 
ence of a Jahn- Teller effect in the cavity QED setting is 
not possible by considering a direct coupling between the 
true atomic ground state and an excited state. A full mi- 
croscopic derivation of the Hamiltonian contains a "self- 
energy" term that has been neglected. This is justified 
in most experimental situations, but in the \g\ > \fuAl/2 
regime it can no longer be left out. Once this term is con- 
sidered, sum-rules or gauge arguments can be employed 
to show that the lower adiabatic potential indeed never 
attains two minima (l7l |. In typical experiments, how- 
ever, the atomic state |<7) is not the true atomic ground 
state. Moreover, by utilizing Raman transitions the cou- 
pling is not direct and the above arguments do not ap- 
ply. Lon ger discussions on this topic can be found in 
Refs. EE 111. 



B. The ex E model 

By considering a degenerate bimodal cavity, the j3 x E 
model can be extended to the e x E one provided the 
atom-field couplings obey d eg ■ e x £ x = hga x and d eg ■ 
e y £y — hg&y. We label the two modes by x and y. The 
Hamiltonian then takes the form 



H eE = tuj\ 



P 2 + P 2 X 2 + Y 2 \ , nn 



2 j-r—a z + hgP x a x +gP y d-y, 

(8) 



with the adiabatic potential surfaces 

(9) 

Instead of having the double well structure as in Fig. [TJ 
the lower adiabatic potential has a sombrero shape. This 
model was extensively studied in Ref. [l6j |. focusing on 



p2 , r>2 [7yz 
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the effects due to the non-zero Berry phase acquired when 
encircling the conical intersection. Initializing a coherent 
state in one of the cavity modes such that it predom- 
inantly populate the lower adiabatic potential, it was 
demonstrated that over longer time periods population 
with be periodically swaped between the two modes and 
that this period greatly depends on the Berry phase. 



C. Renner-Teller model 

If the degeneracy is not conical it is glancing, meaning 
that the tangents of the two potential surfaces are iden- 
tical. This defines the Renner-Teller models. The same 
type of Jahn- Teller effect is possible for Renner-Teller 
models, but the Berry phase vanishes whenever the inter- 
section is in configuration space. From a physical point 
of view, the model we consider is thereby different since 
the intersection occurs in momentum space rather than 
in configurations space. 

To achieve a glancing intersection we look at a three 
level A- atom with the two lower atomic states |1) and 
1 2) coupled via two cavity modes to an excited state |3). 
The model Hamiltonian is taken as 



The general Hamiltonian ([4]) may be rewritten as 



H RT = hw\ . y + 1 



+ 9 (|3)(1|P X + |3)(2|P J/ +H, 



.3=1,2,3 



(10) 

As a three-level system, there are three adiabatic and 
diabatic potential surfaces. Assuming degenerate modes 
and degenerate ground atomic states, E\ = E2 = 0, the 
adiabatic potentials become 



V a ° d (P x ,P y ) = hu;- 



A 2 +P 2 



from which it is clear that V ad (P x ,P y ) and V® d (P x ,P y ) 
possess a glancing intersection at the origin. 



III. EFFECTIVE GAUGE POTENTIALS 

Gauge theories arise in a variety of fields in physics. 
For the most familiar case of a charged particle in an elec- 
tromagnetic field the theory is Abelian. An example of 
a non- Abelian gauge theory is the one of Yang and Mills 
describing strong interaction. Getting an experimental 
handle of a system showing non- Abelian characteristics 
is therefore very attractive and has led to many sugges- 
tions. One possibility is adiabatically evolving systems 
either by means of external changes of the Hamiltonian 
or for ultracold atoms subjected to spatially varying light 
fields 0. 



where 



Pk-A k ) xl 



j 

(12) 



A k = - Eij d ij ■ e k £ k /huj k , $ = -ft Efc UkAl- (13) 

The operators A k and <t have the properties as vector 
and scalar potential respectively. That is, they transform 
appropriately under unitary transformations 20]. For 
any atomic basis, these gauge potentials are matrices and 
they are said to be Abelian if [A k , Ai] = V k and I and 
non-Abelian for non-commuting operators. As we now 
demonstrate, all three examples of the previous section 
render different types of gauge potentials. 



A. Abelian gauge potential 

As the atom-field interaction only includes a single 
mode for the (3 x E model, the vector potential consists 
of a single component which reads 



A = --**. (14) 
Consequently, the gauge potential is Abelian. 



Non-Abelian 5(7(2) gauge potential 



For the second model, the e x E one, we have 



(A x ,A y ) = {a x ,a v ) 

UJ 



(15) 



and hence [^4^,^] = g 2 a z /ui 2 showing that the gauge 
potential is non-Abelian. 

Naturally, the dynamics is considerably richer for a 
system exhibiting non-Abelian structures. In general, 
time-ordering becomes important, e.g. enclosing a "loop" 
clockwise or anti-clockwise will not result in the same 
final system state. In Ref. 20], the time evolution of 
an initial state consisting of one empty mode and the 
other with a coherent state was numerically simulated. 
By properly choosing the phase of the coherent state, 
it will either set off clockwise or anti-clockwise around 
the conical intersection. Figure [2] displays an example 
of the atomic inversion W{t) = (a z ) during such time 
evolutions. The simulation utilizes realistic parameters 
of the qdot cavity QED experimented presented in 0], 
and furthermore takes into account for both cavity losses 
and atomic spontaneous emission. Solid and dotted lines 
corresponds to the different directions around the conical 
intersection of the coherent state. The difference between 
the two curves of the plot demonstrates the non-Abelian 
property. 
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FIG. 2: Time evolved atomic inversion for an initial coher- 
ent state boosted clockwise (dotted) or anti-clockwise (solid) 
around the conical intersection of the e x E model. Parame- 
ters can be found in Ref. [2p| ]. 

C. Non-Abelian SU(3) gauge potential 

The last of our models, the Renner- Teller one, turns 
out to be non-Abelian as well. As a three-level system 
we express the gauge potential in Gell-Mann matrices 

(A x ,A y ) = -°-(\*X), (16) 

UJ 

where A 4 = |3)(1| + |1)(3| and A 6 = |3)(2| + |2)(3|. We 
note that [A^,^] = ig 2 X 2 /uj 2 , with A 2 = -«|1){2| + 
i|2)<l|. 

IV. ANOMALOUS HALL EFFECT 

In the condensed matter community, the coupling of 
our e x E Hamiltonian is said to be on Rashba spin- 
orbit form [2H ]. It is known that such coupling gives 
rise to an effective Lorentz force inducing a transverse 
Hall current. The force is state dependent; opposite sign 
for the appropriate two internal states. When O = 0, 
the population ratio between the two internal states is 
balanced and therefore the Rashba coupling brings about 
a transverse spin current [22j. A non-zero fl breaks this 
symmetry and thereby it causes a net transverse particle 
current. Since the phenomenon derives from an intrinsic 
spin-orbit coupling, and not from an externally applied 
magnetic field as in the regular Hall effect, it has been 
termed anomalous Hall effect [23|. 

The net effective force acting on the state can be de- 
rived from 

dr 1 

F=- = ¥ [H,[H,f}}=g(PxZ)a z , (17) 

where r — (X, Y), P = (P x , P y ) is the momentum oper- 
ators in the Heisenberg representation, and dot indicates 
time derivative. In quantum optics, the field is conve- 
niently examined in phase space. The force l|17j) acts on 
the motion of the phase space distribution. Thereby, the 
Hall current occurs in phase space and hence does not 
involve a true particle current. In the quadrature repre- 
sentation it is clear that any field states evolving under 
the Hamiltonian (|5J| are bounded by a harmonic poten- 
tial. Without the Rashba coupling a Gaussian coherent 



state will bunch back and forth in the two-dimensional 
potential maintaining its chape. On the other hand, if the 
Rashba coupling is non-zero, the motion of the Gaussian 
wave packet will start to bend [24[. Hence, the rocking 
motion will be accompanied by a rotation around the Z- 
axis. This rotation is cither clockwise or anti-clockwise 
depending on the internal atomic state. These assump- 
tions are verified by propagating an initial state with 
mode a; in a coherent state with amplitude a = 10/ v2, 
mode y in vacuum, and the atom in its lower state. The 
results are presented in Fig. [3] showing the trajectories 
of the averages (X) and (Y). It is seen that after hun- 
dreds of oscillations in the harmonic trap, all population 
has been transferred from the x to the y mode. The 
fact that the oscillating amplitude decreases during the 
evolution results from system losses. 
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<X> 

FIG. 3: Time evolution of the expectations (X) and (Y). Due 
to the transverse Hall current, the population initially resid- 
ing in mode x is swaped to the y mode at the final time of 
propagation, which is after approximately 80 ns. The param- 
eters are the same as those used in Fig. [2] and can be found 
in H|. 



V. RELATIVISTIC EFFECTS 

When conical intersections as the one of the e x E 
model appears in momentum space they are frequently 
referred to as Dirac cones. Such Dirac cones have at- 
tracted a great amount of interests in especially the re- 
search of graphene. In the vicinity of the cones, the dis- 
persions are linear as for free relativistic particles. The 
electrons of graphene may therefore show relativistic ef- 
fects like for example Zitterbewegung. For avoided inter- 
sections the relativistic electrons have a non-zero effective 
mass, while for unavoided intersections the electrons are 
massless. 

For small P x and P y we neglect the quadratic terms of 
the e x E Hamiltonian 

Hrei = hg °kPk + h-a z + V{X, Y). (18) 

k— x,y 

This has the form of a Dirac equation for a spin-less parti- 
cle moving within a two-dimensional harmonic potential 
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V{X, Y) [2^1 . The effective masses are equal but with 
opposite signs for the positive and negative energy solu- 
tions. In the present model, small values of P x and P y 
imply that the field amplitudes should be small, some- 
thing easily achieved with ultracold high-Q cavities. 

One interesting observation presented in |26[ is that 
in the non-relativistic limit, the Dirac oscillator in 2+1 
dimensions 



H 



Do 



(p — imfiuur) + [3mc 2 , 



(19) 



with c the speed of light, p and r momentum and po- 
sition respectively, a>k = off — diag(<7fc, &k), and — 
diag(l, — 1), becomes identical to an e x E Jahn- Teller 
Hamiltonian. Somewhat surprising is that the spin-less 
version of Hd , in which the Dirac four-component ma- 
trices afc — ► Ofc and j3 — > <7z, can be mapped onto the 



Jaynes-Cummings model 27] . Furthermore, even in the 
non-relativistic limit the trembling "motion" character- 
izing Zitterbewegung survives and is interpreted as the 
Ramsey interferometric effect (28|. 



ture representation for the fields, the dynamics of the 
combined atom-field system can be thought of as an ar- 
tificial particle moving on a set of coupled potential sur- 
faces. From thereon it is easy to identify different model 
Hamiltonians as Jahn- Teller ones. By rewriting these 
Hamiltonians, we defined effective gauge potentials, both 
Abelian and non-Abelian. Finally, we also demonstrated 
that relativistic effects should appear for weak field am- 
plitudes. 

From the list of references it is clear that it is only 
recently that these systems has been considered as quan- 
tum simulators. A natural conclusion thereby is that 
much more is to be discovered within this topic. To men- 
tion a few possibilities; spin Hall effects in for example the 
Renner-Teller model, Dicke models with multiple number 
of atoms, relativistic effects, spintronics. 

This work did not cover a current proposal for sim- 
ulation of Lipkin-Meshkov-Glick many-body model by 
means of cavity QED [2§| . 



VI. CONCLUSIONS 
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